Diffraction and boundary conditions in semi-classical open billiards 
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The conductance through open quantum dots or quantum billiards shows fluctuations, that can be 
explained as interference between waves following different paths between the leads of the billiard. 
We examine such systems by the use of a semi-classical Green's functions. In this paper we examine 
how the choice of boundary conditions at the lead mouths affect the diffraction. We derive a 
new formula for the S'-matrix element. Finally we compare semi-classical simulations to quantum 
mechanical ones, and show that this new formula yield superior results. 
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Lateral quantum dots, also called quantum billiards 
serve as model systems in the study of the relation be- 
tween quantum and classical physics. A very large num- 
ber of experiments have been donCpOa-.ta'ansport through 
open semiconductor quantum dot, Emma and also exper- 
iments on analog systems such as microwave billiards 
have been performed.ErElCI'B These studies show rapid con- 
ductance oscillations as function of energy or of an ap- 
plied magnetic field. To provide a quantitative descrip- 
tion of these oscillations, several approaches have been 
adopted, such as random matrix th eojy JH nii merical so- 
lution of the Schrodinger equationnE3lI2Hl3lI3 and semi- 
classical(SC) methods. In the SC view, the conductance 
oscillations arise due to interference between pairs of clas- 
sical trajectories, carrying a quantum mechapical(QM) 
phaseJlilOEa between the leads of the billiardB The SC 
approach has been used to describgipStatistical proper- 
ties of the conductance oscillationsMES including weak- 
localization line shapes in chaotic and regular cavities, 
and fractal conductance in systems with mixed phase 
space. Pfie should however be careful in relying on these 
results. E3 The SC approach can also provide an interpre- 
tation of specific frequencies in the conductance oscilla- 
tions, by relating them to specific classical trajectories 
in a billiard. In calculating conductance or transmission 
amplitudes of a system, a SC approximation of the sys- 
tems Green's function is used. In relating this Green's 
function to the transmission amplitude between leads of 
the billiard, one has to take into account diffraction ef- 
fects at the lead mouths. In this paper we study how the 
choice of boundary conditions (BC) on the Green's func- 
tion affect the diffraction and make comparison to QM 
calculations. 

We study a billiard with hard walls and a zero inner 
potential. The Hamiltonian inside the billiard is 
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HW + eA) 2 
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where A is the vector potential of a magnetic field, e 
is the electron charge and m* is the electron's effective 
mass. The dynamics is decided by the Schrodinger equa- 
tion 



(H - E)ip = 0. 



(2) 



We define the Green's function 

(Jf-£)G(q,q')=%,q') 



(3) 



The SC approximation of the Green's function in two 
dimensions istll 



G 6G (q,q',M 



2tt 



^T|A,( q , q ',Ml 1/2 (4) 



(2nih) 3 / 2 



x exp 



-i p (q,q ,k F ) - —\i, v 



where the density of trajectories -D p is taken to vary 
slowly in comparison with the phase from the action S p 
and /i p is called the Maslov index.EJo 

To calculate the S-matrix, we start by letting the 
Hamiltonian act on the wave function ip and expand the 
expression to 
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The divergence theorem states 
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where C is the boundary to the area S, and n' is an 
inward normal to the boundary C. Let 
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We also get 
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Equations 
tity 



dq" 



results in an analog to Green's first iden- 
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We take the complex conjugate of equation ([Uj), inter- 
changing <p and tp, and subtracting it from equation (10) 
after some manipulations we get 
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We would however like to get rid of the magnetic terms 
on the left hand side, we now insert 



eA 

F = — i>(P* 
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(12) 



into the divergence theorem and choose a gauge such that 
V • A = 0, resulting in 



dq' — hpA • V</>* + <p*A • W] 
s m 
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dl'—A-ti^*. (13) 
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This allows us to rewrite equation (|ll[), while replacing 
0*(q / ) = G(q,q') as 
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a version of Green's theorem. 

We now calculate the diffractiop-,from a lead mouth, 
following in the steps of Kirchhoff.E 2 ] The area S is taken 
to be the half circle, see figure 0a. The contour G is taken 
to be divided into two parts Ci, the entrance lead mouth 
and Co, the rest of C. We make the requirement that 
tp = and dtp/dn' = along the wall, this requirement 
is elaborated on later in the text. With zero magnetic 
field, G, ip ~ r^ 1 / 2 , and the integration along the half 
circle part of C tends to zero as the radius r — > oo.E2l 
The motivation for letting r — > oo within a finite billiard 
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FIG. 1: a) Geometry of the entrance lead, S is the half circle 
enclosed by C = Co U C\ , where Ci across the entrance only, 
n' is an inward normal to C, a' is an angle at the entrance, q 
is a coordinate and x and y' form a local coordinate system, 
b) Geometry of the exit lead, n is into the lead, x and y form 
a local coordinate system. Path C2 is across the exit and C3 
is along the entire wall. 



is that in the SC theory the wave function along all clas- 
sical trajectories is the free space wave function. The 
total wave function is then the superposition of all such 
trajectory wave functions. In nonzero magnetic field we 
are dealing with edge states, localized near the wall. If 
we assume that the current is going in the — y-direction, 
then ^>(q') will vanish everywhere on the half circle ex- 
cept near the wall in the lower half-plane (y' < 0). On 
the other hand the Green's function G(q, q') is also lo- 
calized to near the wall but for G(y ^> y',y'), ie. the 
Green's function at the interior point q will not see an 
excitation at a point q' that is far enough downstream 
the wall. For an interior point q either ip(q') or G(q, q') 
will go (exponentially) fast to zero on the half circle when 
r — > 00. It is thus enough to integrate over C\ instead of 
G in equation (|l~4|). Using equations @ and (||) in jl^ ) 
we can obtain the wave function inside the billiard 
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where dy' = —dl' . 

To calculate transmission amplitudes, wer-sjtart from 
the definition of the current density operatorcil 



J(<l) 



2m* 



VS(q-q') + S(q-q')V 



(16) 



where V = —ihV + eA is the momentum operator. We 
obtain an operator for the current into the exit lead by 
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integrating across it 



J = [ dyj{y)-h, 
Jc 2 



(17) 



where Ci is across the exit and y the transverse coordi- 
nate, see figure |l]b. The eigenmodes n in the entrance 
lead are taken to be 



<(*'.!/) = 



/\ Ak n x 



where v n is the velocity in the lead, k n the k-vector, and 
x' and y' are defined in figure The eigenmode wave 
functions of the exit ip^x^y) are defined analogously. 
The continuation ip n of ip„ into the billiard is given by 
equation (|l^). The scattering matrix element is equal to 
the matrix element 
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By inserting equations (|l5|) and ([l8]) we get 
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+ ^ (mv m + 2eA(y) • x) — 

+-o (™ m + 2eA(y) • x) (mv n + 2eA(y') ■ x') G 

(20) 

where primed coordinates relate to the entrance and un- 
primed to the exit. 

Equation ( |20| ) is the transmission amplitude calculated 
using the Kirchhoff approximation, which to recall con- 
sists of the following assumptions: 

1. ip and dtp/On vanish everywhere on C except on 
the lead mouth. 

2. The values of ip and dip/dn on the lead mouth are 
equal to the values of the incident wave in the lead. 

For example Schwieters et al.0 use this approximation. 
There are however some inconsistencies in the Kirchhoff 
approximation.E3 It can be shown for the Schrddinger 
equation that if ip = and dip/dn' = on any finite 
surface, the the only solution is ip — everywhere. This 
inconsistency can be lifted by choosing Green's functions 
with appropriate BC, i.e. either Dirichlet or Neumann. 
We start with the Neumann Green's functions 



^(q,q')=0forq'on C. 



(21) 



If the contour C is an infinite line across the entrance and 
S all space to the right, the Neumann Green's function 
can be obtained by the method of images 



Gjv(q,q') = G(q,q')+G(q,q") 



(22) 



where q" is the mirror image of q', reflected in C. In the 
case where q' lies on C, the Neumann Green's function 
reduces to 

G jV (q,q')=2G(q,q'). (23) 
We then get the wave function in the billiard to be 



iP = - I dl'G N (q,q') 



h 2 dip 
2m* dn' 



ip A ■ n 

m 



(24) 

The coupling to the exit lead should be handled in a way, 
analog to the coupling to the entrance lead. The equation 
for the scattering matrix element should be symmetric 
with respect to the direction, renaming the entrance and 
exit leads should make no difference. We therefore need 
to make the same restriction as on the entrance lead, 
i.e. we make a demand on the Green's function on a 
boundary G3 crossing the exit 



dG 



N 



dn 



(q,q') = for q on G 3 



This is again done by the method of images 

G jV2 (q,q') = G jV (q,q') + G J v(q" , ,q') 



(25) 



(26) 



where q'" is a mirror image of q, reflected in G3. For 
q'eC and q G G3, we get 



G jY2 (q,q')=4G(q,q'). 



(27) 



Using this Green's function and equation (g4j) in equation 
(O), we obtain 



m y/V m V„ 



dy / dy'C m {y)Uy') 



(mv m + 2eA{y) ■ x) {mv n + 2eA(y') ■ x') G, (28) 

which in case of zero magnetic field or if gauge can be 
chosen such that A • n = and A • n' = 0, reduces to 



ihyJVmVn f dy I dy'Cn(y)^n(y')G, (29) 



which we recognize as the QM expression for the S- 
matrix. This expression is also most commonly used for 
SC calculations. However in QM calculations the Green's 
function is derived in a way that includes the effects of the 
leads. This expression is exact when using a QM Green's 
function. This is not true for the SC Green's function 
which totally ignores the eigenmodes of the leads. Given 
that the wave function ip must be zero on the billiard 
walls, the choice of Neumann BC is obviously wrong, 
since G partially coincides with the billiard wall. We 
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conclude that this commonly used expression does not 
use the correct BC. 

The correct BC to use with the SC Green's function is 
Dirichlet, 



Gn(q,q') = for q' on C. 



(30) 



The Dirichlet Green's function can also be obtained by 
the method of images 



G D (q,q')=G(q,q')-G(q,q") 



(31) 



The normal derivative of the Green's function on C can 
be computed to 



dG D (q,q') =2 _dG(q,q') 



(32) 



dn' dn' ' 

for q' on C. The wave function in the billiard is in this 

case 



2m* 



/ lQ dG D (q,q r ) 



dy i>' 



Ci 



dx' 



(33) 



Following the path of derivation for the Neumann Green's 
function we arrive at 



q 



<h, I dyT m (y)Uv')^ (34) 



which is our new formula for calculating S'-matrix ele- 
ments. 

Details on calculating trajectories, for the SC Green's 
function can be found in references We need however 
to take the derivative of the Green's function too. We 
will here restrict our self to zero or weak magnetic field, 
because the derivatives of the semi-classical Green's func- 
tion in magnetic field are very complicated to compute. 
In zero magnetic field the action only depends on the 
length of the trajectory and it is enough to only look at 
the ends of the trajectory to calculate a derivative of the 
action. In non-zero magnetic field the action will depend 
on the total geometry of the trajectory and thus also the 
derivative. In the following derivation we assume zero 
magnetic field, and in the end of this article we show that 
the resulting formula still yield very good correspondence 
to QM, even in finite magnetic field. The density of tra- 
jectories D p in equation (Q) is considered to vary slowly 
in comparison to the phase given by the exponential func- 
tion. The contribution of one trajectory to the Green's 
function is considered to be a plane wave in the vicinity 
of the trajectory. The <9/<9a;'-derivative then reduces to 
taking the derivative along the trajectory and projecting 
on the x'-direction, ie. multiplying with sina^, where 
the angle a' s of the trajectory is defined in figure |l|a. 



dG SC (q,q',k F ) 
dx' 

h~~~ s dq< 



2tt 



^2 lA^q'^ 



,1/2 



i . l S p (q,q',k F ) 
x — sina — - — — exp 



p 

-b p {q,q,k F ) - —fi p 



(35) 



where q'^ is a coordinate along the trajectory. The action 
can be written as 



Sp — Tik F l P j 



(36) 



where l p is the length of the trajectory. We can then 
calculate a derivative as 



dSp 



-Tiki 



(37) 



and the Green's function derivative 



dG SC (q,q',k F ) 



dx' 



(27ri^) 3 /2 

sin a' s exp 



1 1/2 



-6 p (q,q ,k F ) - 



(38) 



The above expression for the Green's function derivative 
and the corresponding derivative of x have been used 
in the following transmission computations. We have, 
as mentioned above, used equation ( |38"| ) and the corre- 
sponding ^-derivative as approximations in weak mag- 
netic fields, we have although inserted the exact expres- 
sion for the action S p into these derivative approxima- 
tions. This results in an exact phase but approximate 
amplitude for the contribution of a trajectory. The re- 
sulting good agreement with QM, as will be seen below, 
justifies this approximation. 

Calculations were made on more than ten different bil- 
liards (triangles and rectangles with different aspect ratio 
and different lead positions) of which two are included in 
this article, see insets in figures || and The compu- 
tations have been made both using QM for comparison 
and using SC theory with the three different equations 
(|2"o|), (|2S|) and ([hL. Details on QM calculations can be 
found in rcfcrcncco. The conductance through the two 
billiards have been analyzed in terms of their length spec- 
trum of the transmission/reflection amplitude defined as 



t(l) = dk F t(k F ,B = hk F /{er c ))e~ 



ilk f. 



(39) 



where r c = hk F /(eB) is the cyclotron radius, and t is the 
transmission amplitude, ie. an S-matrix element. The 
length spectrum is thus in magnetic field calculated at 
constant cyclotron radius, so the geometry of the trajec- 
tories stays constant. The action along a trajectory in a 
magnetic field can be written as 



S p (q, q', k F ) = hk F (l p + A p r c 1 ), 



(40) 



where A p is an area related to the trajectory. With the 
density of trajectories varying slowly with respect to k F 
we can see that the length spectrum of the transmission 
amplitude will be strongly peaked at I = l p + A v r~ x , or in 
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FIG. 2: The transmission and reflection length spectrum for Sn calculated QM and SC using Kirchhoff approximation [eq. 
(pfj|)1, Neumann BC [eq. (0)] and Dirichlet BC [eq. ([m])] for a rectangle, see inset, in zero magnetic field 
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FIG. 3: The transmission length spectrum for Su calculated QM and SC using Kirchhoff approximation [eq. 
BC [eq. (p8|)] and Dirichlet BC [eq. (p^)] for a triangle, see inset, in a magnetic field with cyclotron radius r c 
is side of triangle. 
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zero magnetic field I = l p which makes the name "length 
spectrum" obvious. 

The two analyzed billiards included in this paper, see 
insets in figures || and ||, are a rectangle, with sides of 
length L and 1.1L in zero magnetic field and a trian- 
gle in magnetic field with r c = 0.6L, where L is the 
side of the triangle. The leads have a width w = L/15, 
and the transmission and reflection amplitudes has been 
calculated in the regime of one conducting mode in the 
leads. There is a good correspondence between the posi- 
tions of the peaks in QM and SC calculations. The peaks 
that are missing in SC calculations are due to "ghost" - 
trajectories, which bounce diffractively against the lead 
mouths-JIhcse trajectories are not included in this SC 
model.E3n2l The included billiard geometries have been 
specially chosen to avoid that "ghost" -trajectories coin- 
cide in length with normal trajectories which otherwise 
makes comparison with QM computations difficult. 

There is a very good correspondence between SC cal- 
culations using Dirichlet BC [eq. (p4[) l and the QM calcu- 
lations. The Neumann BC [eq. Q2S|)j still gets the peaks 
in the right positions but can not reproduce the correct 
amplitude. The Kirchhoff approximation [eq. (p0|)] can 
be seen as a mean value between Neumann and Dirichlet 
BC and produce peak amplitudes between the Neumann 
and Dirichlet amplitudes, ft is interesting to note that 
even in the case of the triangle where the geometry near 
one of the leads differs from the case in the derivation 
and we are using an only approximate expression for the 
Green's function derivative, the Dirichlet BC still gives 
superior results (there is no derivative in the expression 
for the Neumann BC transmission amplitude, so it is not 
affected by the approximation). The actual difference 
between the usage of the three different BC is, at least 
in zero magnetic field, in the diffraction from the leads. 
The different BC yield different diffraction patterns, see 
figure U 

To conclude, we have studied the effect of the choice 
of boundary conditions at the lead mouths. We find that 
in contrast to prior practise of using an expression from 
QM that agrees with Neumann BC, or using the Kirch- 



hoff approximation, one should use Dirichlet BC. This is 
not in disagreement with QM because the QM Green's 
function differs from its SC counterpart in that it includes 
the effect of the eigenmodes in the leads, while the SC 
Green's function does not. The SC Green's function is 
only an approximation of the QM Green's function and 
because of its ignorance of the wave function outside of 
the billiard it should be treated differently from the QM 
Green's function. We therefore propose a new expression 
for the scattering matrix, equation (|3~4]). By compari- 
son to QM calculations, we show this expression to yield 

superior results. 

r 1 1 TT7- 1 1 r~t 




FIG. 4: The Fraunhofer diffraction pattern from a single lead 
for kF = 1.5ir/w, where w is the width of the lead. The 
plot shows the amplitude of the wave functions in the far- 
field as function of entrance angle a, with the three different 
boundary conditions. 
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